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Abstract 

Possible strong Van der Waals force is searched in the nuclear interaction, whose 
existence is expected when the fundamental interaction between quarks is strong or 
superstrong Coulombic type. The relation between the type of the extra singularity 
at t = of the scattering amplitude and the power of the asymptotic behavior of the 
long range potential is derived. The square root cusp in the once subtracted S-wave 
amplitude, whose existence is expected when the Van der Waals force of the London 
type is acting, is observed when we analyze the high precision phase shift data of the 
low energy proton-proton scattering. 
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1 Van der Waals force in the composite model of hadron 



Before we accepted the composite model of hadron in 1960's, the nucleons were the 
elementary particles and the nuclear force was essentially the short range force arising 
from the exchanges of a pion and a set of pions and heavier particles. However the 
short range nature of the nuclear force is not necessarily correct in composite hadron, 
especially when the fundamental force which combines the quarks is Coulombic. It is 
strange that even after the acceptance of the composite model, the nuclear force is still 
regarded as short range with argument that, when the momentum transfer is not very 
large and one does not explore the inside of the composite system, we can regard the 
nucleon as elementary particle with good approximation. However this argument is not 
true, especially when the fundamental force is Coulombic. In such a case, as we shall 
see later, the Van der Waals interaction between hadrons must always appear. It is 
instructive to review the mechanism of the appearance of the Van der Waals potential 
between the neutral 'atoms' and to estimate the strength C of the potential in the 
asymptotic region 



where o is the radius of the composite particle. 

Let us consider a simple system in which two atoms are in the ground states of 
L = 0, and they are fixed with separation R. Because of the quantum fluctuation, 
atoms can also stay in the excited L = 1 state, whose duration is At ~ h/AEi by the 
uncertainty principle where AEi is the excitation energy. Suppose that atom 1 is in 
the excited P-state and therefore has the dipole moment d\. The dipole must produce 
the field G at the location of atom 2, which decreases as R~ s for large R. The field G 
induces a dipole d 2 of the atom 2, which is a 2 G where a 2 is the polar izability of the 
atom 2. The induced dipole d 2 produces in turn a field G' at the location of atom 1, 
and therefore the energy of the system changes by —d\ ■ G' which decreases as R~ e for 
large separation R. This is the mechanism of the appearance of the attractive Van der 
Waals potential of the London type. In particular for the case of the ordinary atoms, 
since the electric charge is the source of the constructive Coulomb field, the dipole d 
is the electric dipole and the field G is the electric field. However the mechanism to 
induce the Van der Waals force is more general, and any charge, which produce the 
fundamental Coulombic force, can do the same job. For example, the color charge in 
QCD or the magnetic charge in the dyon modelpQ is responsible for the appearance of 
the Van der Waals interaction between hadrons. 

In order to estimate the strength C of the Van der Waals potential V(R) ~ —C/R®, 
let us consider the simplest case in which atom 1 and atom 2 are respectively composed 
of two particles of opposite sign. As shown in fig.l, fi is the relative coordinate of atom 
i. The interaction Hamiltonian H 1 of the system of atom 1 and atom 2 is 
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Fig.l. System of atom 1 (a- A) and 
atom 2 (b-B) separated by R. 
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, where *e 2 is the "fine structure constant" of the relevant charge which produces the 
fundamental Coulombic field. The energy shift V(R) of the two-atom system is 

If we see that the numerator of each term of the summation of Eq.(3) is positive 
definite and the energy denominators are not less than the first excitation energy AE± = 
Ei—Eq, we can obtain the lower and the upper bounds of the coefficient C of the Van der 
Waals potential by retaining only the first term and by replacing all the denominators 
by AEi respectively, and which is 

\<l\(R^H')\0>\^ c ^ <0\(R^Hr\0> ) (4) 



AEi ~ ~ A£i 

in the estimation of the upper bound the closure relation is used. If we substitute 
Eq.(l-2), the upper and the lower bounds become 

-r\r^ < C < -— -ri T2 , (5) 



3 (AEx) 1 2 ~ ~ 3 (AEi) 

where fj 2 is the ordinary mean square radius of the ground state of atom i. On the 
other hand, r is the radius, which relates to the transition amplitude, defined by 

f 2 = 3 1 5| 2 with z= I d 3 ripl i ^rcosO^ofifi . (6) 



Now we can estimate the order of magnitude of the strength C of the Van der Waals 
potential, in which AE\ = 2 and the radius is 1/2 ~ 1/6 in the unit of the pion mass 
or the pion Compton wave length. The largest value of the radius 1/2 comes from 
the charge radius of the nucleon, whereas the smallest radius 1/6 is obtained from the 
distance of the repulsive core 1/3 of the nucleon- nucleon interaction. It is interesting 
that in QCD, in which *e 2 ~ 0.3, the strength C becomes 0.003 even if the radius 
is 1/2. On the other hand in the magnetic monopole model of hadron,[Q in which 
*e 2 ~ 137.0/4, C becomes order of 1 when the radius is l/4~l/5, and in this case the 
Van der Waals interaction can compete with the nuclear force arising from the meson 
exchanges. 
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2 Can the nuclear potential be derived from the meson 
theory ? 

In 1950's and 60's, the nuclear potential had been studied expensively, in which Take- 
tani's idea played an important roll. In his idea the nuclear force is divided into three 
parts according to the range. The outermost region is described by the one-pion ex- 
change, whereas the innermost region is treated by phenomenological way and the hard 
core potential is used in \iR < x c . On the other hand in the middle region, the two-pion 
exchange is responsible for the nuclear potential, and which is supposed to be calcu- 
lated from the interaction Lagrangian L % n of the meson theory. In figure 2 the central 
potential 1 V r ( t(i?) of the singlet even state is plotted against R. The curve OPE is the 
one-pion exchange contribution. The curve TMO is the OPE plus the two-pion ex- 
change contribution which is computed by the standard perturbation (static), whereas 
HM is the same quantity with full recoil considered. |2j 




(b) (a) 

Figure 2: Potential 1 Vq(R) is plotted against R in the unit of the neutral pion mass and its Compton 
wave length. The curves a(0.30), b(0.25), c(0.20) and d(0 ~ 0.15) are the one-boson exchange potentials 
of Eq.(2-1) fitted to the experimental data with various hard core radii, which are given in the brackets. 
The curve OPE is the one-pion exchange potential. The curve TMO is the OPE plus the two-pion 
exchange potential which is computed by the perturbation theory (static), whereas the curve HM is 
the same type potential with full recoils are involved. Some attractive potential seems to be missing 
from the theoretical curves (TMO and HM) in the region around R ~ 1. 
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On the other hand, curves a(0.30), b(0.25), c(0.20) and d(0 ~ 0.15) are the one- 
boson exchange potentials: 

| oo for a: < x c 

in which a; = uR and the parameters yl and B are determined by fitting to the phase 
shift data, and whose values are listed in the table below. The four curves differ by the 
hard core radii x c , and whose values are written in the bracket. 



curve 


x c 


A 


B 


(a) 


0.30 


-12.1 


16.2 


(b) 


0.25 


-13.6 


23.0 


(c) 


0.20 


-14.2 


26.1 


(d) 


~ 0.15 


-14.8 


28.2 



It is evident that large attractive potential is missing from the TMO and HM 
potentials, and therefore these theoretical potentials can not reproduce the phase shift. 
In order to settle the discrepancy, it is customary to add the term of the fictitious 
<7-meson exchange, in which the a-meson mass m a and the <r-N coupling constant are 
free parameters. However in this paper we shall examine whether we can attribute 
the missing attractive force to the strong Van der Waals interaction mentioned in the 
previous section. 



3 Extra singularity of the amplitude characteristic to the 
long range interaction 

In order to identify the long range interaction in the nuclear force, it is important to 
find the behavior of the scattering amplitude A(s, t) characteristic to the existence of 
the long range interaction. When the nuclear force arises from the meson exchanges, 
A(s,t) has a pole at t = 1 and the continuous spectrum At(s,t) starting from t = 4 
and so on. What is important is that in such a case A(s,t) is regular at t = 0. On the 
other hand, when the long range force is acting, an extra singularity appears at t = 
and the spectral function becomes non-zero from t = 0. Since t = is the end point 
of the physical region — Au < t < 0, where experimental data are available, we can in 
priciple determine the type and the magnitude of the singularity at t = 0, when the 
precise data are available in the small neighborhood of t = 0. The threshold behavior 
of the spectral function A t is specified by the power 7 and the coefficient C, which are 
defined by 

A t (s,t) = irC'f(l + o(t)) . (8) 

These parameters relate to the parameters a and C of the potential V(R) which are 
defined by V(R) ~ —C/R a in the asymptotic region. The relations are 

2C' 

a = 2 7 + 3 and C = -^-T{2- 1 + 2) , (9) 
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where m is the reduced mass of the potential scattering. [3] 

It is not difficult to derive the relations, if we remember that the potential V{R) 
can be represented as the superposition of the Yukawa potentials: 

V(R) = ~ / w(t')e~ RV¥ dt' . (10) 

7rm z H Jo 

Since the range of the integration is from to infinity, this is nothing but the Laplace 
transformation of RV{R). What is important is that the weight function w(t) is the 
imaginary part of the first Born amplitude (s,t). we can see this simply by comput- 
ing the Fourier transformation of both sides of Eq.(10) and which leads to a dispersion 
relation of the form 

TT Jo V — t 

In particular if the weight function is w(t') = 7rCY 7 , the potential of Eq.(10) becomes 
or" i 

y(i?) = -^ r(27 + 2) i?W • ( 12 ) 

Since the change of the weight function at finite t' does not alter the tail of the potential, 
we obtain the relations between the powers a and 7 and also between the coefficients 
C and C given in Eq.(9). 

Therefore the extra singularity of the scattering amplitude A(s,t) at t = is 

A{s,t) = -^—C'{-tr + - , (13) 
sin 7T7 

where dots mean the background regular function and can be represented by a polyno- 
mial of t. When 7 is an integer re, we must take the limit 7 — > n, and the singular term 
becomes C't n \og{—t). Since t = —2v{\ — z), there are two ways to observe the extra 
singularity at t = 0. One is to see the angular distribution for fixed v and to observe 
the singular behavior (1 — z) 7 at z = 1. Second is to make the partial wave projection 
and to observe the singularity v" 1 at the threshold in the partial wave amplitude ai{v). 
In our normalization, a^v) relate to the phase shifts in the elastic region by 

a e (u) = ^=e^ , (14) 

or in terms of the effective range function X^(i/) 

aAv) = —j= where XAv) = — ^cotcWz;) . (15) 

*/(")- ^* 

In order to make the extra singularity more visible in the search we shall use the once 
subtracted partial wave amplitudes which are defined by 

ar(i/) = «0M-"0(0) and a once {v) = Mv± f Qr ^ 
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In ag nce (v), the singularity is changed to C"z/ Y_1 . Therefore the extra singularity of 
the Van der Waals interaction of the London type must appear as the square root cusp 
C'^v 1 ! 2 at v = . Moreover since the Van der Waals force is attractive, namely C > 0, 
the cusp must point downward. 

In order to observe the extra singularity as clearly as possible, it is necessary to 
remove the nearby singularities and prepare the domain of analyticity of the back 
ground function as wide as possible. First of all we must remove the unitarity cut, 
whose spectral function is Im af (v) = ( v / i/2y / I7) sin 2 S^u), or in terms of the effective 
range function X^(u) 

to «« w = ^WSA with *M-^«**<"> • < 17 > 

This is done by introducing a function 

Kr e (») = aT ce (v) - - fdv' -^ Q f , (18) 

vr Jo vV(z/ - v) 

in which the domain of the integration [0, U] may be chosen as the region where the 
reasonably accurate phase shift data are available. 

The amplitude K^{y) is sometimes called Kantor amplitude, which is free from the 
singularity at v = when all the forces are short range. Concerning the left hand cuts, 
the nearest one is the cut of the one-pion exchange which starts at v = (—1/4). Such 
a cut is removed in (K$ nce (v) - a° £ nce ' l7T »), where 

«r'"M = (i£)i(^(i+^)-*.«) ■ (") 

The next nearest left hand singularity appears at v = — 1 which is the threshold of 
the two-pion exchange. Since the threshold behavior of the spectrum is 1 — 

v) 3 / 2 9(— 1 — u), the distortion of the background function arising from the two-pion 
exchange in the small neighborhood of v = is small. Therefore we can observe the 
extra singularity at v = unless its coefficient is not very small. 



4 Observation of the cusp at v = in the once subtracted 
S-wave amplitude 

In order to distinguish the effect of the long range force from that of the two-pion ex- 
change, we need accurate data. It is well known that the accuracy of the measurements 
of the low energy proton-proton scattering is highest in the hadron physics. Therefore 
we shall analyze the once subtracted S-wave amplitude of the p-p scattering aQ nce (v), 
which is spin-singlet. Because of the electromagnetic interaction between protons, some 
modifications are necessary when we construct the once sutracted Kantor amplitude 
KQ nce (u) given in Eq.(18). Contrary to the phase shift 8q{v) of the scattering by purely 
short range potential, in which the incoming and the outgoing waves are the partial 
wave projection of the plane wave, we must use the electromagnetically distorted wave 
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as the initial and the final states. The information of the scattering of £ = state is 
summarized in the phase shift Sq(u), 5q(u) or Sq (v). The small differences of these 
phase shifts are due to what portion of the electromagnetic interaction is involved when 
we compute the distorted incoming and outgoing waves. The potentials which distort 
the initial and the final waves are: (C) pure Coulomb force of the point proton, (E) 
the Coulomb force of the proton with electric form factor plus the vacuum polarization 
effect or (EM) in addition to the interactions of the case E other tiny effects of the 
electromagnetic interactions such as the interaction between the anomalous magnetic 
moments of protons are involved, respectively. 

The relation between the effective range function Xq{v) and the phase shift of 
Eq.(15) must be modified. Its explicit form is given in other paper. What is 
important is that the effective range function Xq(v) is regular at v = 0. If we notice 
this property, it is straightforward to construct the once subtracted Kantor amplitude 
KQ nce (v) which is free from the singularity at u = when the interactions are short 
range plus electromagnetic. [3] In order to ease the search of the extra singularity, it is 
necessary to subtract the term of the one-pion exchange at least approximately. The 
S-wave phase shift data 5q M {u) of the energy dependent analysis by Nijmegen groupjS] 
are used in the evaluation of the once subtracted Kantor amplitude Ko nce (u). In figure 



-Kl nce (v) = -(Kg nce (u) 



onceAir , 
l 



v)) is plotted against I} a b, in which the range of 



Ti a b is < Ti a b < 125MeV.. Figure 4 is the same graph with enlarged scale and the 
narrow energy range < < 20MeV.. In figure 4, the error bars are attached to 
the data points (circles) of very low energy region, however for higher energy the error 
becomes smaller than the size of the data point. We can clearly see a cusp at v = 
pointing upward, which means the attractive long range force is acting. 




20 40 60 80 TU0t9£5 



T(MeV.) 




Figure 3: ~Ko nce (u) is plotted against T iah in 
Ti a b < 125 MeV., in which the Nijmegen phase 
shift data of the energy dependent analysis are 
used. The curve is the fit by the spectrum of the 
long range force in the range 0.6 MeV. < Ti a b < 
125 MeV.. The graph indicates the square root 
cusp at Ti a b = 0, which points upward. 



Figure 4: The enlarged graph of -KS nce (v) is 
plotted in T ia b < 20 MeV. The three curves 
are 3-parameter fits to the data in 0.6 MeV. < 
Tab < 125 MeV. The curve (L) 3 is the fit by 
the spectrum of the long range force, whereas 
the curves of (s 0)3 (dots) and (s 6)3 (dash) are 
the fits by the short range potentials. 
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In order to determine the type of the long range force, let us make the chi-square 
fit by the spectrum of the long range force 



A t (s,t) =nC'f<e- pt 



(20) 



Best fit in 0.6MeV. < T lab < 125MeV. is achieved by 



7 = 1.54 



C = 0.18 and (3 = 0.063 



(21) 



in the unit of the neutral pion mass. The curve of the best fit is shown in fig. 3 and also 
(L)3 curve of fig. 4, which designates the long range spectrum with three parameters. 
On the other hand, the fits by the short range force, whose spectrum is the sum of 
three delta functions and their coefficients are free parameters, are shown in fig. 4. For 
the dot curve (50)3, the locations of the delta functions are t = 4, 9 and 16, whereas 
for the dash curve (sb)^, the locations are t = 9, 16 and 25 respectively. It is evident 
that the short range force cannot reproduce the data points well. 

5 Conclusion and prospect 

Since the parameters of the spectrum of the long range force are determined in Eq.(21), 
by using the relation Eq.(9) we can evaluate the parameters of the tail of the long range 
the potential V{R) ~ —C/R a , which become a = 6.09 and C = 0.18 in the unit of 
the neutral pion mass. It indicates that reasonably strong Van der Waals force of the 
London type is acting between nucleons. Therefore contrary to the ordinary model of 
the nuclear potential, in which only the short range forces are involved, it is important 
to construct nuclear potential anew based on the strong Van der Waals interaction plus 
the short range potential arising from the exchanges of mesons. 

If we remember the fact that the Van der Waals interaction is universal, the extra 
singularity at t = must occur in every amplitudes of the hadron-hadron reactions. 
We can expect to observe it whenever very accurate data are available, or when we can 
construct a back ground function with wider domain of analyticity. The low energy 
p-p amplitude is the former case, whereas the low energy P-wave amplitude of the 
pion-pion scattering belongs to the latter case. Because of the crossing symmetry we 
can remove the cut of the two-pion exchange as well as the unitarity cut from the tt-tt 
amplitude, therefore we can expect to observe the cusp in a\{y)/v when reasonably 
accurate pion-pion phase shifts data are given. Details are found in separate papers. |Hj 

Finally from the upper bound of the strength C of the Van der Waals potential 
given in Eq.(7), we can estimate the lower bound of the strength *e 2 of the fundamental 
Coulomb interaction, it becomes 



in which the values of the nucleon radius f = 0.5 and of the first excitation energy 
AEi = niA — ni-N = 2.17 are used as well as the strength C = 0.18 of the Van der 
Waals potential. It is interesting to compare the result Eq.(22) with the strength of the 
QCD, which is *e 2 ~ 0.3, and also with that of the magnetic monopole model of hadron, 



*e 2 > 3.06 



(22) 
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which is *e 2 = 137.0/4 by Dirac's charge quantization condition. Therefore the value 
of the strength C of the Van der Waals potntial supports the magnetic monopole model 
of hadron rather than the QCD. My dream is to confirm the magnetic momopole model 
of hadron 1 by observing directly the monopoles of opposite sign fuse to form a meson. 
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